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Abstract 

We analyze multi-instanton sector in two dimensional U(iV) Yang-Mills the- 
ory on a sphere. We obtain a contour intregral representation of the multi- 
instanton amplitude and find "neutral" configurations of the even number instan- 
tons are dominant in the large N limit. Using this representation, we calculate 
the 1,2,3,4 bodies interactions and the free energies for N = 3,4,5 numerically 
and find that in fact the multi-instanton interaction effect essentially contribute 
to the large N phase transition discovered by Douglas and Kazakov. 



1 Introduction 



In the last two years, there is remarkable progress toward understanding relation be- 
tween large N QCD and string theory. Important contribution to this subject was 
made by Gross and Taylor [Q]. They showed that "two dimensional QCD is a string 
theory" . Strictly speaking, the partition function of U(iV) Yang-Mills theory on Rie- 
mann surface E^ is equivalent to that of no fold string theory with target space Eg, in 
other words, sum over branched covering maps of Riemann surface E^. But Douglas 
and Kazakov discovered that on the topology of sphere the system has a 3rd order 
phase transition at (XA) C = tx 2 in the large N limit and showed the equivalence is 
restricted within the strong coupling region XA > ix 2 0. 

Originally the phase transition was discovered from the analysis of so called Cauchy 
problem in this system. It is well known that when we calculate matrix model partition 
function in the large iV limit, we encounter the Cauchy problem in eigenvalue distribu- 
tion function. In our case the partition function has the similar structure as the one of 
Gaussian matrix model and we encounter a two-cut Cauchy problem in A A > tx 2 . By 
exploring the free energy, there is 3rd order gap between the weak and strong coupling 
region. 

One physical explanation of the phase transition is the following. From the string 
point of view (it is view from one side i.e. strong coupling region), covering map from 
sphere to sphere (which corresponde tol/N leading term of Yang-Mills free energy on 
sphere) allows arbitrary but even number of branch points (2n-2) according to arbitrary 
number of winding (n). Each branch point has competing entropy (target space area 
A) and energy. In additon there are many ways to construct world sheet from n 
target spaces. If the entropy beat the energy , the string expansion becomes infinite 
and the phase transition occurs ||. This scenario succeed to explain the large N phase 
transition qualitatively and semi-quantitatively(it gives critical coupling (XA) C ~ 11.9). 
But this point of view seems to be insufficient because it lacks global perspective and 
gauge field point of view. 

Recentry, Gross and Matytsin showed that instanton induces the phase transition 
0. They obtained the following results. 

1. Yang-Mills partition function on sphere can be expressed by sum over instantons. 

2. The instanton sector gives the weak coupling result. 

3. In the weak coupling region the 1 instanton amplitude of charge ±1 has exponen- 
tial damping factor and thus is negligigle. 

In addition, the 1 instanton amplitude of charge ±1 becomes oscilating in the strong 
coupling region. But it's contribution to the free energy are order N~ 1 / 2 and thus 
negligible in the large TV limit. Hence it is important to investigate the effect of the 
multi-instanton. 
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In this paper we explore the effect of multi-instanton and find out that there is 
some sort of neutrality in the large N limit. That is, "neutral" configurations of the 
multi-instanton are dominated in this limit. We also find that the interaction rather 
than the entropy essentially cause the phase transition. We conjecture 2 bodies force, 
4 bodies force, 6 bodies force,. . . equally contribute to the large N phase transition. 

The content of this paper is following. In section 2 the multi-instanton amplitude 
and the partition function are rewritten to multiple contour integrals and the neutrality 
is found. In section 3 nature of one instanton and two instantons systems are analized in 
detail. In section 4 numerical calculation of the multi-instanton amplitude and the free 
energy are performed. In appendix using the above formulation, the multi-instanton 
contribution to Wilson loop is rewritten to multiple contour integrals. 



2 Contour integral representation and neutrality 

We begin with Migdal-Rusakov heat kernel representation for the Yang-Mills partition 
function on a sphere @, 

r N t 

Zqcd = J 2M M exp(- — J d 2 x^tiF^F^) 

= Y:(dimi?) 2 exp(-Mc 2 ( J R)), (1) 

where A is area of the sphere, R is irreducible representation of the gauge group 
U(iV) and C 2 (R) is the value of the second Casimir operator of rep R. Because R 
is characterlized by Young tableau, sum over row lengths in Young tableau replaces 
sum over representations. We remark that the condition that Young tableau must be 
stairway -like is irrevant because dimi? has form of the Van der Monde determinant 
and thus has permutation symmetry. Then the sum becomes free sum over Z N and we 
can rewrite it by the Poisson resumation formula. We get the following expression [J] 
ill: 

Zqcd = e^^^f 2 £ e^^^M), (2) 

N 



II d Vi 11(4 - 4vr 2 mJ)e"^ 

1=1 1<J 

= / Hdy i A(y i + 2nm i )A(y l -2nm i )e-^^ (3) 

J — OC ■ i 



1 Hereafter we absorb A into A 

2 A similar expression was obtained by M. Caselle et al ||. They showed that the phase transition 



is due to the winding modes of "fermion on circle" . 
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where y^ = yi~yj and = mi — rrij. In ||, |], [5| it was showed that {m} corresponde 
to all Euclidean classical solutions (which we call instanton) up to gauge transforma- 
tion and has nonperturbative effect with respect to both 1/N and A. Hereafter we call 
mi instanton charge and call number of nonzero m^s instanton number. The w({m}) 
looks like Gaussian matrix model but there is a deformation in the Van der Monde de- 
terminant A. In this section we rewrite w({m}) using the method of ortho-polynomial 
(in this case Hermite polynomial) and obtain a new multiple contour integrals repre- 
sentation. The new representation makes clear the role of the multi-instanton which 
was not known. 



2.1 and 1 instanton sector 

In the case of instanton (i.e all = 0),w({m}) is itself parition function of Gaussian 
matrix model and has the following form: 



"(0) = j_J{dy^\y)e-^^=Z G {—)=c N {—y^. (4) 



In the case of 1 instanton (i.e. rrii = m5ik), Using the property of the Van der 
Monde determinant and the Taylor series expansion of the Hermite polynomial, we 
obtain, 

r°° JL __N_ y^iV 2 

w(m) = / [[dyiA(yi + 2irm5i k )A(yi - 2irm5 ik )e ^ 2^=i^ 
J -°°i=i 

N ^ r°° _ al 2 

= (JV-l)! 22 h ---h n ---h N -i / dy k P n {y k + 2nm)P n {y k - 2vrm)e 
n=o J -°° 

i N-1 n 



1 „ , N s r dt i^m 2 Nt , 1, 



iV ^ A n =0 «=0 ft 



where 



and 



N GK 2A J J 2ni v t> ' v 7 



AT 

P ™( x ) = 7^7WTW H n(\l ~ x ) ( 6 ) 



/°° ^P n (y)P m (y)e-^^ 2 = Mnm = v / 2T(-) n+ ^!5 nm . (7) 

This contour integral representation for the 1 instanton amplitude was obtained in || 
from another aspect. They calculated it using large N saddle point method and found 
it to have behavior change from exponential damping to oscilating at A = m 2 Ti 2 . The 
detail is given in section 3. 
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2.2 multi-instanton sector 

In the same way, we can rewrite multi-instanton amplitudes to multiple contour inte- 
grals: 

w({m}) 

| [ dyiAfa + 2mm i )A(y i - 2nmi)e ™ 2^=i Vi 
■°°i=i 

f°° N N _JV V-iV 2 

yi s § n ^ s / n^n p ^w(^+ 27rm i) p ^«(^~ 27rm i) e ^^=1^,(8) 



fj,eS N o-eSjv 00 i=i i=i 



where Sjv is the permutation group on iV object. 

According to number of nonzero m's , i.e. number of instantons, the following 
analysis is different. Assume number of nonzero m's is k, using permutation symmetry, 
nonzero m's are driven to m/s from i=l to i=k. Then we have, 



N 

w({m}) = S SW £ II h a (i)8<T(i),„oa(i) 

/iG5jv ctSSjv i=k+l 



k f°° N 2 

x II / dyi p v(i){yi + 2irm i )P fl0(T ( i) (yi - 27rm i )e~3z 1 '* 

i=1 J-oo 

= (N - k)\ J2 sg 11 ^ S h ---h ai ---h ak ---h N _ 1 

k 

t 00 N 2 

x I] / dyiPaXVi + 27rm i )P lx(ai) (y i - 2irnii)e~^ . (9) 

i=1 J-oo 

In the same way as 1 instanton sector, using the Taylor series expansion of the Hermite 
polynomial, we obtain for the above multiple integrals, 

k roo N 2 

1 1 / dViPaiiVi + 2nm i )P fl{ar) {y i - 2-Km^%~^^ x (h ai ■ ■ ■ h^)' 1 

J -co 



1=1 

k min(a ( ,//(a;)) /m O-rrmAT 

= TT V (f!^il)a,-i, ( _ 27rm) Ma,)-i». (10) 

The each series in the above equation can be expressed as a contour integral by the 
following transformation formula. 



min(a,b) 
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where the contour of t encircles the origin counterclockwise. Using this formula, we 
obtain w({m}) for k instantons, 

= N(N - 1) • • • (N ~k + 1) „,|i_§; gn "/ Si ' "t 2^1— t 

)<c -^(m?t 1+ -+m?A)( m Ml) 1 +*Ml) )oi . . . ( "VW L±^W )"fc.(12) 

We remark the configurations such as a\ = a 2 do not affect the above equation. Hence 
we can drive J2a^---^a k to independent sum ' ' 111 that expression there 

are many pole free terms, then , after eliminating pole free terms, this integrals finally 
become 



w({m}) 



Z G (ij) i dh r dt k 



N(N - 1) • • • (N - k + 1) J 2ni J 2m 
x e -^("^i+---+»^ fc )(i + . . . (i + L)" 

The last sum is just determinant of matrix M which has (ij) element 

Ma = 777-- . (14) 

It's diagonal element is 1 and this determinant has all information about the interaction 
between instantons. 



2.3 partition function 

We can rewrite the partition function as another contour integral. For this purpose, we 
must expand w({m}) by number of instantons and classify elements of the permutation 
group Sfc into the conjugacy classes. By symmetry argument, the contribution from 
different elements belonging to same conjugacy class is same. Then we obtain, 

z QCD = z w ™ k j2 E e -^ (m ? + --- +m ' ) / — ••• /— e- JV( * m i ( * i)+ --- + * m " ( *" )) 

n=i) mi,...,m n fO 

x Y, sgn[a]T[a}M lai M 2a2 ■ ■ ■ M nan , (15) 



conj. class 



where 



Z wcak = e^- 1 )(^) Jva Z G (^)=oonste ,va ^-i k « A )-«, (16) 
^ . . 4ir 2 m 2 t , .„ 1, 

S m (f) = _-__log(l + -). (17) 
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Z weak is the partiton function in the weak coupling region and it has simple A depen- 
dence. We assume a has cycle structure [1 CT1 • • •n <Tn ](a\ + • • • + na n = n) and define 
T[a] as number of elements in the conjugacy class which a belongs to, 

T[a] 



sgna = (_)^+<x4+- ) ( 19 ) 

E = E • (20) 

conjclass <ri ,...,<r n — 



J"XH Ynan— n 



We use (12345) (67) type element as the representative element of the conjugacy class. 
Because the multiple integrals factorize according to the cycle structure of a, contribu- 
tion from the multi-instanton are exponentiated with a constraint (o"i + • • • + na n = n) 
which is expressed by another contour integral. We get the following result: 

Zqcd = z~ f £ exp(£ tp^), (21) 



3=1 



where 



J ^ ,„ J 2m J 2m 

m\,...,m,j^p\) 



xM 12 M 23 ---M jl (22) 
is the "connected" amplitude of the j instantons. 



Some remark should be noted. In the large N limit, the multiple contour integrals 
for w({m}) or <x,- are dominated by the saddle points. Fortunately, it is no need for 
solving complicated saddle point equations in our case. The saddle point equations 
are decoupled and the solutions of the equations are copies of the saddle point found 
by Gross and Matytsin. Then if we neglect the interaction term M12M23 • • • Mji, atj 
becomes factorized such as (aiY and the behavior change from exponential damping 
to oscilating at A = TT 2 m 2 observed in || still hold for the multi-instanton amplitude. 
But there is a remarkable feature in the interaction term of the multi-instanton sector. 

We point out that in the region of A > n 2 special charge configurations of the 
multi-instanton exist. If we consider the following configurations in ctj (even j): 

mi = —ni2 = = —7774 — ■ ■ ■ — m, (23) 
the large N saddle points for any t{ are same and satisfy the next relation: 
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When A > m 2 TT 2 , both saddle points must be selected. Moreover when tj = t± and 
U+i — tzp, all Mjj + i's have pole and usual N~ 1 / 2 power law (Gauss integral factor) must 
be changed. Toy example of such phenomena is the following: 

l-expf-AW))^, (25) 

where t c is saddle point of 4>(t). We can approximate this to a principal value integaral 
along the steepest descent line and a contour integral along small semi circle around 
the saddle point and obtain finite value with order exp(— N(f)(t c )). We call this 
type of problem "singular saddle point". Same power change occurs in the multi- 
instanton sector. The configuration that all Mu + i have pole is dominant contribution 
with respect to N in the strong coupling region. We call this phenomena "large N 
neutrality" . By "neutrality" , we would not intend general neutral configuration which 
satisfy X^Li m i — 0, but intend that all instantons have equal absolute charge and 
number of positive charge instantons and number of negative charge instantons are 
same. Hence by definition "neutrality" in aj with odd j is ruled out. Bellow we will 
see aj with odd j is in fact negligible in the large N limit. 

3 One body and two bodies effective action 

Below we will see generel 2n bodies interactions equally contribute to the phase tran- 
sition. Even so, it would be worthwhile to analyze familiar two bodies interaction 
between instantons in detail. We find that the interaction between instantons is sensi- 
tive to A and m. 

3.1 one body effective action 

First we consider the contribution of one instanton with charge m to the partition 
function. In the large N limit, the saddle point equation $' m (t) = has two solution 

t± = l± ^ 1 m , We see quite different behavior according to A^7r 2 m 2 . 
3.1.1 m 2 > A/n 2 

Both saddle points exist on negative real axis and satisfy 

$ m (f_) < $ m (t+). ( 26 ) 

Neverthless only t + saddle point is selected since t_ steepest decent line is not consistent 
with the original contour. Hence we obtain in the large iV limit , 

w(m) ~ (-) 1 ^= e - JV * m(t+) . (27) 

v /2vriV|$^ ) (t + )| 
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If we define j(x) as || 

x , 1 + >/T _ •'' 



7 (x) = v/T^ - - log (28) 

2 1 — Vl — ^ 



we obtain 



o 2 2„ I \JV-1 



Z linst = e — tw(m) = ; 1 J e — (29) 

Since 7 (a;) is positive real for x < 1 but pure imaginary for x > 1 , self energy term 
survives in this region but disappears at A = m 2 TT 2 . 

3.1.2 m 2 < A/n 2 

In this case both saddle points are selected because t± steepst descent lines are con- 
sistent with original contour and satisfy 

Re[$ m (t+)] = Re[$ m (t_)]. (30) 

Therefore, we get 

Z lmst = ( ~ )JV = (e %-^^> + e-** ^^-r(^)). (31) 

Since j(x) are pure imaginary in this region, the above equation has order N' 1 ^ 2 in 
the large N limit. 

3.2 two bodies effective action 

Next we consider the contribution of two instantons with chage m and m! to the 
partition function in the large N limit. 

3.2.1 m 2 > m' 2 > A/n 2 

Straight forward application of large N saddle point method gives 

1 1 



e - JV (*«.(^)+* m '(»+))detM(t + , s + ), (32) 



, , -l±Jl-A/Tr 2 m 2 -l±Jl-A/TT 2 m' 2 , 

where t± = — 5 , s± = — 5 and 



1 

detM(t, s) = 1 - — — — — ; -. (33) 
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As a further approximation we assume m' 2 ^> A/tt 2 , then we obtain 

^ ]_ e -2N(q2+q' 2 )+N(logq 2 +\ogq' 2 )+\og(q-q>) 2 ^ ^ 

where q = ^n 2 /Am is effective charge and we used the following expansion formula of 
7 around x = 0: 



+ ilog(i) _ i+s ^_ (i 



) n+1 . (35) 



We note that the interaction term is next leading order in 1/N and the self energy 
term strongly suppresses this configuration. 

3.2.2 m 2 < m' 2 < A/tt 2 and m ^ —mf 
In the same way we obtain, 

w(m,m') = 1 = 1 = e i f"- jV( ' I> - (t + )+ ^ (s + )) detM(t + ,g + ) 

y/2TrN\$£\t + )\,j2nN\Q%(s + )\ 

+ 1 1 e-^*^-> + ^ s +»detM(t_, s + ) 

y/2irN\^(t.)\yj2irN\^(s + )\ 

+ -, 1 = . 1 =e- JV($m(t+)+ ^ ,(a - )) detM(t + ,s_) (36) 
^ttTVI^^I^ttTVI^^)! 

+ 1 . 1 e-'l 7r - Jv ^^-) + ^^-»detM(t_, s_). 

^ttTVI^^)! V2ttJV|$S(s-)I 

If we further assume m' 2 <C then we obtain 

^2inst ~ ^[0(qq') + ^(-gg')e log{9 - <? ' )2 - los{( ' +9 ' )2 ], (37) 
where we used the following property of 7 : 

7(x) ~ -—x x > 1. (38) 

We note that in this region the self energy term almost disappears and another inter- 
action term log (9 + q') 2 is added. 
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3.2.3 m 2 < A/tt 2 and to = -to' 

As is alreasdy noted, in this case we encounter the singular saddle point. 



-to 



dt_ r ds_ N(» m( * )+ * m , W ) (1 _ 1 y (3g) 

2mJ 2m K (1 + t + s) 2 ' y ' 



Since the 1st term of the determinant has no singularity at saddle points, we can neglect 
this term and obtain 



w 



(to, -to) ~ -N 2 r dzze- Nz f{z) 2 , (40) 
Jo 



where 



/(z)E /^ e_( ^ +,)J¥t(i+ ! )ff - (4l) 

From above equation f(z) is essntially ol\ with minor changes. If the z is above 4 — 
47r 2 m 2 /A , exp(— Nz/2)f(z) has the behavior of exponential damp. Hence we can 
restict the z integration range up to 4 — An 2 m 2 /A. In addition, exp(— Nz)f(z) 2 has 
both oscilating and constant terms . Oscilating term corresponde to t — s — t± 
which give order N^ 1 for exp(— Aw ™ N )w(m, —to), and constant term correspondeto 
t = t±, s = t T which give order for exp ( — 4n2 ™ 2 N ) w (m, —to). In the large iV limit 
we obtain 

2N r it ,2ir 2 m 2 . 2n 2 m 2 l~A~ . 

^2inst--— x - arcsin — — - 1) —J——-1 . (42) 

tt 2 A A V vr z m z 

4 Numerical calculation 

In this section we show the result of numerical calculation of a±, a2, a^, and free 
energies for A^=3,4,5 using the Mathematica . We assume < A < An 2 . From the 
above argument, in this region we can restrict the instanton charge as rrii = ±1 (which 
we call 'truncated') in multi-instanton sector. 

Figl,2 show that iV = 10, a = A/n 2 = 2 results respectively for ai,a 2 ,a 3 and a 4 . 
We observe that compared to the a 's with a > 1 , the a 's with a < 1 are negligible. 
Moreover we see remarkable difference between a±, (even j) and a 2 , «4 (odd j) with 
respect to order and shape in the region of a > 1. The shapes of a a are almost same 
and both have the properties of non-oscilating and linear scaling with respect to N. 
On the other hand the shapes of a±, a% are different. But both have the properties of 
oscilating and N~~ l l 2 scaling . We already observed that ot\ has N~ 1 ^ 2 scaling and a 2 
has linear(iV) scaling using the large iV saddle point method. The numerical calculation 
indicates the agreement with eq (29), (42). From this observation and the argument of 
section 2, we conjecture the following scenario of the large iV phase transition. That is, 
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otj with even j have linear scaling with respect to iV in the region of a > 1 and inside 
the exponential in eq (21), we have 

N [N/2] 1 

E K - A —j* 3 E * 0(N 2 ). (43) 

j=l J n=l Z71 

This gives order N 2 result for the free energy in a > 1. Thus aj with even j equally 
contrubute to the large N phase transition. 
Fig3 show the truncated free energies 

1 i ri dz 1 1 ,^ i m ^ 

= iV 2 " l0g[ / 2™ ^Tl^ 6Xp( g ^ 

for iV=3,4 and 5. The partiton function for N=5 turns to have negative value in a^3.2. 
This causes the singular behavior around a ~ 3.2 of the N = 5 graph (here we plot real 
part of the F). We observe that the graphs stand up at a ~ 1 and in the neighborhood 
of a = 1 the graphs seem to converge in the large iV limit. But in the region a^2.5 
discrepancy of the iV=3,4 and 5 graphs is large. We think the truncation becomes 
invalid in this region. Hence we may conclude that near A = it 2 nature of the large 
N phase transition remain at iV = 3, 4, 5 level, but we cannot distiguish whether this 
behavior change is smoose or not. 

5 Outlook 

We have reformulated the multi-instanton amplitude as the multiple contour integrals 
and found that this represenation make clear the properties and the effect of multi- 
instanton in two dimensional QCD. In particular "neutral" configurations of the even 
number instantons are dominated in the large N limit. But we think further investi- 
gation is needed about 

1. derivation of the order of large N phase transition in terms of our contour integral 
representation, 

2. clearcut order parameter and Ginzburg-Landau effective action, 

3. possibility of phase transitions at A = n 2 k 2 (k :positive integer). 
Some explanations are needed. 

About point 1, we did not obtain the free energy analytically because we do not 
know general method for solving the singular saddle point problem. But if such method 
is found out, we will obtain analytic form of the free energy and hence the order of the 
phase transition. 

The reader may be confused with resspect to point 2. By 'clearcut', we intend order 
parameter which is zero or nonzero according to A^n 2 . In order to drive argument 
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of phase transition to Ginzburg -Landau type, we must know such order parameter a 
priori as variable of effective action. In our case instanton number is naive contender, 
but it is ruled out. From the expression of instanton amplitudes they are not positive 
definit and can not have the meaning of probability^ ( Boltzman weight). Hence we 
must consider other order parameter. If such order parameter is found out, the effective 
action will be the following form near critical point 0. 

r(0) = -(/5-/3 c )0 2 + ^ 3 (g>0) (45) 

About point 3, we ignore the instanton configuration of m — ±2, ±3, ... in sec- 
tion 4 since such configuration give exponentioal damping factor in the region of 
7r 2 < A < Att 2 . Therefore, same argument leads to that in the region of k 2 7r 2 < 
A < (k + l) 2 7r 2 (k : positive integer) we can ignore the instanton configurations of 
m = ±(k + 1), ±(k + 2), . . .. Hence the multi-instanton configurations which should be 
included are changed at A = k 2 ir 2 . Other approaches did not predict phase transitions 
at A = k 2 ir 2 , but we think whether the changes lead to some phase transitions or not 
is still important problem. 

We hope this work will shed some lights on the global structure of two dimensional 
QCD. 
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3 That is why we don't use the word 'instanton condensation'. 

4 In the case of n th order phase transition, g<j} 3 term is replaced by gcj) 2n ^' 
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Appendix Multi-instanton effect in Wilson loop 

In this appendix we note the multi-instanton effect in Wilson loop. Gross and 
Matytsin also showed the formula for Wilson loop on a sphere [[/J: 

W n (A u A 2 ) = {±trU n ). (46) 

Here the loop devide the sphere with area A into two disks with area Ai, A 2 and U is 
gauge field holonomy along the loop . 

W n (A 1 ,A 2 )Z QCD 

a' n 



' V k=l {m} 

inAo „ s . , _ inA\ 



x 



f°° T-r , a / mAo r~ \ a / „ vnAi r . iv ,,2 

/ J] d vMVi + + —±5 ik )A( yi - 2<n mi + — _L<y tt ) e -az 2^=i ^ . (47) 

J— oo • -, iv iv 



From above equation, calculating Wilson loop is almost same as calculating the parti- 
tion function with additional imaginary charged instanton. Hence in the same way as 
the partition function, we can rewrite this as multiple contour integrals: 

= ^wcaky^I v e _^ (m 2 + ... +m?) rdti..(dh 
f-i.ll ^ , n J 2m J 2m 

xe -^Kti+-+mf t! )(l + 1)JV . . . ^ + I)iV 
r± e -2ninrn 1 ^+2^nrn 1 ^^t 1 -4^t ldet (l) B 



_L /^±i e -^^-(l + -Vdet^C], (4* 
JM J 2m U + i 



where 



and 



Ba = -3 , (49) 

^ = (50) 



A = -2irm i + —A 1 6 il , (51) 
7i = —2mrii (for i = l, ...,n), (52) 

7n+i = jjA x . (53) 
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Using this formula we can calculate multi- instanton effect of the Wilson loop in the 
large N limit. 

For example in instanton sector (weak coupling phase), we obtain, 

W n (A u A 2 )Z QCD \ 0inst = Z weak i f ^e-^(l + j) N . (54) 
By rescaling t — > Nt, we obtain the known result: 



W n (A 1 ,A 2 ) ~ i,/-i-J 1 (2nJ^). (55) 



n V A^ 2 V ^4 
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Figure captions 

Fig 1: The truncated a 's for N = 10 and various a = A/ir 2 are shown. 

Fig 2: The truncated a 's for a = 2 and various N are shown. 

Fig 3: The truncated free energies for iV = 3, 4, 5 and various a are shown. 
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Fig3:The truncated Free Energies for N=3,4,5 
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Fig1 : The truncated ALPHA'S for N=10 
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Fig2:The truncated ALPHA'S for a=2 
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